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Abstract 

We review and discuss some different techniques for describing 
local dispersion properties in fluids. A recent Lagrangian diagnostics, 
based on the Finite Scale Lyapunov Exponent (FSLE), is presented 
and compared to the Finite Time Lyapunov Exponent (FTLE), and 
to the Okubo- Weiss (OW) and Hua-Klein (HK) criteria. We show 
that the OW and HK are a limiting case of the FTLE, and that 
the FSLE is the most efficient method for detecting the presence of 
cross-stream barriers. We illustrate our findings by considering two 
examples of geophysical interest: a kinematic meandering jet model, 
and Lagrangian tracers advected by stratospheric circulation. 
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1 Introduction 



Transport processes play a crucial role in many natural phenomena. Among 
the many examples, we just mention particle transport in geophysical flows 
which is of great interest for atmospheric and oceanic issues. The most nat- 
ural framework for investigating such processes is the Lagrangian viewpoint, 
in which the tracer trajectory x(t) is advected by a given Eulerian velocity 
field u(x, t) according to the differential equation 

f = u(x.«) (1) 

Despite of the simple formal relation (|ID, the problem of connecting the Eu- 
lerian properties of u(x, t) to the Lagrangian properties of the trajectories 
x(t), and viceversa, is a very difficult task. In the last 20 — 30 years the 
scenario has become even more complex by the recognition of the ubiquity 
of Lagrangian chaos (chaotic advection). Even very simple Eulerian fields 
can generate unpredictable Lagrangian trajectories which are practically in- 
distinguishable from those obtained in a complex, turbulent, flow |^, ^, |]. 

In the following we will restrict our attention to the case of two-dimensional 
incompressible velocity field u(x, t), with x = {x,y). Incompressibility is au- 
tomatically satisfied by introducing the stream function = ?/;(x, t) and 
consequently defining the velocity field in terms of partial derivatives as 
u = {ipy, —ipx)- The evolution equation (|1|) then becomes 

dx dy 

Formally (|]) is a Hamiltonian system with Hamiltonian ?/'(x, t). Chaotic 
trajectories x(t) typically appear as a consequence of the time dependence 
of^ i- 

Many geophysical flows, when observed at sufficiently large scale, are 
within this class [^, ^. Moreover time dependence can be often considered a 
perturbation over a given steady flow, i.e. @ represents a quasi-integrable 
Hamiltonian system. It is well known that in quasi-integrable Hamiltonian 
systems chaos can be quite non uniform in the phase space, due to the pres- 
ence of the invariant KAM tori with a chaotic layer around them The 
presence of these regular islands (also called coherent structures in the con- 
text of geophysical flows) is of particular importance for the dispersion pro- 
cess because they act as barriers to transport. The sensitivity of different 
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diagnostics of transport to the presence of barriers will be the main topic of 
our investigation. In particular, we will consider the Finite Time Lyapunov 
Exponent (FTLE), the Okubo- Weiss (OW) and Hua-Klein criteria (HK), and 
the local Finite Scale Lyapunov Exponent (FSLE). 

We discuss these different methods by considering two examples. First, 
we study transport properties in a kinematic meandering jet model, formerly 
introduced for describing the Gulf stream |]^, Second, we analyze a large 
number of stratospheric isoentropic trajectories, computed according to (j^) 
from assimilated wind fields, describing Lagrangian motion around the polar 
vortex. In both situations our results show that the existence of barriers 
limiting the dispersion across the stream is well described by the FSLE but 
it is completely missed by the OK and HK criteria, which can only depict the 
landscape of alternating unstable hyperbolic and stable elliptic points of the 
flow. The FTLE will be discussed in relation to the OW and HK exponents 
and to the FSLE. 

The remaining of this paper is organized as follows. In Section |^ we 
introduce and discuss the different diagnostics for characterizing dispersion. 
Section ^ is devoted to the analysis of the meandering jet and of the polar 
vortex. Finally, in Section ^ we present some conclusions. 



2 Characterization of local transport proper- 
ties 

In presence of Lagrangian chaos, two close trajectories typically separate ex- 
ponentially in time 0]. Thus the natural statistics we adopt for describing 
chaotic particle spreading is the relative dispersion statistics. Relative sep- 
aration between two particles R(t) = x (t) — x(t) evolves according to the 
velocity difference 

— = u(x(t) + R(t), t) - u(x(t), t) . (3) 

As far as particle separation remains much smaller than the typical length 
scale Ie of the velocity field, we can linearize (|^) around the trajectory x 
and, for a generic time dependent flow, we expect exponential growth of the 
separation, i.e. 

R{t) ~ i?(0)e^* (4) 
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where A is the Lagrangian Lyapunov Exponent (LE). 

In the opposite hmit, Ie, the two particles feel uncorrelated velocity 
fields and one recovers the standard diffusive regime, i.e. 

~ 2Dt (5) 

where the average is taken over many particle pairs and where D is the 
diffusion coefficient. 

It is important to remark that, in most realistic situations, both asymp- 
totic regimes, i.e. very small R{t) for (^) and very large R{t) for (^) cannot 
be attained. From one side, particles separation can be not sufficiently small 
to justify the linearization leading to (^. In the opposite limit, large sep- 
arations cannot be reached in presence of boundaries at scales comparable 
with Ie- As a consequence, the asymptotic quantities as A and D cannot be 
computed and a non-asymptotic characterization of transport is needed [§]. 

Let us discuss, now, some techniques one can use to characterize local 
dynamical properties of a system, in particular the relative dispersion rate 
as a function of the initial position. 

2.1 The Finite Time Lyapunov Exponent 

Let us start from the definition of the Lagrangian Lyapunov Exponent: 

A = lim lim - In \ (6) 

t^ooi?(o)^o t R{0) ^ ' 

In d^) one basically assumes that the linearization of the perturbation Rif) on 
a generic reference trajectory holds for an infinite time. This is correct only 
if the perturbation can be considered infinitesimal at any time. The char- 
acteristic time naturally associated to the LE is known as the predictability 
time T\ = A~^, which is the characteristic time at which one can predict the 
position of the tracer in the future. The FTLE is obtained by avoiding the 
limit t — CX3 in (|^). This gives the instantaneous growth rate over a finite 
interval r as 

I, R{t + T) 

7^ t = lim -In^— 7 

^ R{t)^0 T R{t) ^ ' 

which, at variance with A, depends on the initial point x(t). The FTLE is 
distributed around a mean value which is nothing but the LE, < 7,- >= A, 
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where the average is computed over a virtually infinite number of r intervals 
along the trajectory [|ll[]. 

In principle, even at very small R{t), one must wait a certain time interval, 



T^, needed to drive the perturbation along the Lyapunov direction [|T2|. This 
is necessary if we want to measure intrinsic properties of the system. In 
presence of many degrees of freedom, the possibility that the waiting time 

is of the same order of the predictability time Tx cannot be excluded (see 
[y] for a review on the predictability time in extended systems). 

Let us now discuss some practical shortcomings arising when we want to 
analyze realistic situations described by experimental or model data. First, 
as we already rescaled, in quasi-integrable Hamiltonian systems, different re- 
gions in the phase space can display different behaviors . As a consequence 
one has a non trivial spatial distribution of Lyapunov exponents: zero if the 
trajectory lies in a regular island, positive if it diffuses across the stochas- 
tic layer. In special cases, when for instance structures of the velocity field 
are "localized" persistent features, at least within time intervals considerably 
longer than the characteristic Lagrangian time, a more refined description in 
terms of finite-time Lyapunov exponent can be more appropriate. 

In order to measure 7r(^) at a point x(t) one can make use of the fol- 
lowing procedure. Backward integration in time for an interval T* bring the 
trajectory at the point x(t — T^,). An infinitesimal perturbation 5x(t — T^,) is 
switched on and it is integrated forward to (5x(t). It T* is sufficiently long, i.e. 

> Tyj, the perturbation 5x(t) will be aligned along the Lyapunov eigen- 
vector and a further integration to 5x(t + r) will give the FTLE according 



to |]T2[. In general, is not known a priori and can vary with the initial 



conditions. More serious problems arise from the limits of resolution in the 
knowledge of experimental or simulated Lagrangian data, which can easily 
disrupt the linear approximation scenario. 

When R{t) attains finite sizes, i.e. is of the order of the characteristic 
lengths of the system, the so-called Finite-Scale Lyapunov Exponent (FSLE) 
give an appropriate description of the intrinsic physical properties of disper- 
sion at different scales of motion. We will discuss this point in Section 

2.2 The Okubo- Weiss and Hua-Klein criteria 

In two-dimensional turbulent flows, the stirring properties of initially small 
material lines are related to the combined effect of eddy and jet features of 
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the velocity field. In cases when continuous velocity fields are given, a pop- 
ular way used to characterize the local rate of separation of initially close 
trajectories is the Okubo- Weiss (OW) criterion [III, based on the com- 



putation of the eigenvalues of the velocity gradient tensor. If explicit time 



dependence cannot be neglected, Hua and Klein (HK) |]TOl have proposed a 
generalization of the OW criterion, based on the computation of the eigen- 
value of the acceleration gradient tensor, related to the distribution of the 
pressure field. 

Let us recall the two criteria in the case of 2D incompressible velocity field 
with Lagrangian evolution given by (^. The evolution of an infinitesimal 
separation R(t) is given in the tangent space as 



dt V -Ipxx -tpxy 



R^MR, (8) 



where the Jacobian M has the property = AqI with Aq = —det(WL). At 
small times, the solution of (BD is 



R(t) 



1 + Mt + -m^ 



R(0) + O(t') (9) 



2 

where N = AqI -I- dNl/dt. The Okubo- Weiss criterion consists in computing 
Ao, i.e. the product of the eigenvalues$$eigenvalues of M. We recall that the 
quantity Aq can be written in terms of the square strain and the square 
vorticity a;^ as 

\, = -y-u') (10) 

If Ao is positive, the two eigenvalues of M are real, the velocity field is 
locally hyperbolic and strain overcomes rotation. For negative Aq, we have 
imaginary eigenvalues and a predominance of rotation over strain. 

Of course, the Okubo- Weiss criterion may not be sufficient to determine 
the local strain- vorticity balance in a time dependent fiow. In this respect, 
the Hua-Klein criterion, being based on the sign of the largest eigenvalue of 
the N matrix A+ = Aq + Ai, with 




d^XX dipyy 

dt dt 



(11) 



gives a "more Lagrangian" description. Both of them provide a picture of 
the distribution of stable elliptic points and unstable hyperbolic points in the 
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flow. Let us observe that in the case of stationary velocity field, dlSA/dt = 
and one has A+ = Aq. 

Let us remark the relationship existing between the HK criterion and the 
FTLE. The "instantaneous" Lyapunov exponent 7,- estimates the growth 
rate of a typical perturbation within a finite time interval r, after the pertur- 
bation has aligned along the most unstable direction, so that it is an intrinsic 
property shared by all the trajectories (except for a set of zero probability 
measure). The HK eigenvalue A+ estimates the local maximum strain rate, 
regardless any alignment time of the perturbation, so that, from the La- 
grangian point of view, it may be biased by transient behaviors. In other 
words, measuring A+ corresponds to measuring 7,- at very small r (the in- 
tegration time step) starting with a perturbation always aligned along the 
locally most unstable direction. In practical situations it can happen that 
the HK eigenvalue and the FTLE give similar local descriptions but, from a 
theoretical point of view, starting with a perturbation along the local most 
unstable direction is as arbitrary as choosing any other direction: after a 
transient, the time evolution will drive the (infinitesimal) perturbation defi- 
nitely along the most unstable Lyapunov eigenvector. 

2.3 The local Finite Scale Lyapunov Exponent 

In most cases of interest, the linear regime during which the exponential 
growth of the inter-particle distance occurs can be not resolved. Particle 
spreading is generally observed on large spatial scales, of the order of the 
characteristic lengths of the system, and appropriate non linear techniques 
must be employed to quantify relative dispersion rates (see [0 for a review 
about non asymptotic properties of transport and mixing). 

Let us consider a very small (infinitesimal) initial perturbation -R(O) on a 
trajectory x. For a chaotic system, after the initial transient, R{t) typically 
grows exponentially in time according to (^). We wait until R{t) reaches a 
certain threshold at a given time tj. Let Xj the position of the trajectory 
X at time ti. At a later time tf, R{t) will reach a larger threshold 6/ = r ■ Si 
with an assigned r > 1. We now define the r— amplification time of R{t), 
relatively to the initial position Xj, as T(5j, Xj, r) = tf — ti. From this quantity. 
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the local Finite Scale Lyapunov Exponent is defined as: 

Xr{5i,^i) = —^ ^Inr (12) 

The exponent Xr{Si, Xj) is a measure of the local amplification rate of a per- 
turbation of size Si on a trajectory passing by the point Xj. In the case of 
periodic time dependence in the equation of motion, the local FSLE will 
depend also on the initial phase. 

The above described prescription is necessary in order to characterize the 
perturbation growth as an intrinsic property of the system. Of course, in 
realistic situations, we have to consider the following problems: the fiow is 
not simply periodic and thus may depend explicitly on time. Moreover, 
we cannot observe distances below a certain threshold because of finite reso- 
lution, and thus the initial perturbation cannot be considered infinitesimal. 
Another important remark is that in practical applications it may not be 
possible to have information on a uniform distribution of initial conditions. 
In particular, let us consider a 2-D time dependent velocity field (e.g. the 
surface circulation of a sea or an isoentropic layer of the stratosphere) and 
the relative Lagrangian dispersion of trajectories. Usually one cannot set 
the initial distance between two particles to an arbitrarily small size, and 
wait until the relative position vector is aligned with the most unstable di- 
rection and starts expanding with a typical rate measured by the LE. Finite 
resolution imposes a lower limit to physically reasonable distance sizes, say 
5, and we can only hope to take into account all the possible realizations 
of the local dispersion rate by taking the average of Ar(5, x) (r > 1) over a 
large number of directions around the initial point x, i.e. on a sphere with 
radius 6. Furthermore, if 6 is not very small if compared to the characteristic 
lengths of the system, the linear regime of instability is already expired. The 
FSLE, by its nature, is a non linear indicator of trajectory instability so it 
can measure relative dispersion rates at finite scales of motion. How long 
can we follow two trajectories so that their FSLE is still meaningful as a 
local diagnostics ? The answer depends essentially on how much rapidly the 
velocity field varies in time relatively to the Lagrangian characteristic time. 
For instance, if Tg is the time scale within which the Eulerian structures, e.g. 
current systems, change their geometrical and physical aspects, local FSLE's 
are meaningful only if they are observed on times <C Te, i.e. in an almost 
frozen field approximation. 
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Let us now briefly discuss the main objective of tlie present paper, i.e. tlie 
detection of barriers in particle transport. We seed the flow with an uniform 
distribution of Lagrangian tracers and compute the FSLE for any trajectory 
according to the prescription given above. 

If £e represents some characteristic length in the flow, e.g. the typical 
size of the eddies around the Gulf Stream in the North Atlantic Ocean or 
around the polar jet current in the stratosphere, then particles trapped inside 
vortices or traveling down jet streams may never separates beyond the scale 
Ie, giving \r{5) = for 5 < On the contrary, particles located in the 
chaotic layer will give a positive FSLE. Thus at an appropriate value of the 
threshold 5, the map of FSLE ([l^) can be used as an efficient indicator of 
the presence of barriers in transport. 

We will see how a simple periodically perturbed meandering jet can be a 
significant test to show that a more appropriate diagnostics is required when 
macroscopic barriers are under investigation. An interesting result related to 
geophysical data analysis concerns the barrier effects of the jet current of the 
stratospheric polar vortex (southern hemisphere). This technique is being 
used also to study local mixing properties in ocean systems. 



3 Results 

3.1 Numerical experiments 

We first discuss a simple, but not trivial, kinematic model in which a barrier 
to motion is known to exist for certain values of the parameters. The tran- 
sition to the barrier-breakdown occurs by variation of some parameters. We 
use this model to show what kind of information can be extracted from the 
different techniques previously discussed. The system, formerly introduced 
as a model of transport across the Gulf Stream ||, consists a time periodic 
streamline pattern forming an oscillating meandering (westerly) current with 
recirculations along its boundaries: 



\& = —tank 



y — Bcoskx 



+ cy (13) 



where k is the spatial wave number of the structure, c is the retrograde 
velocity in the "far field" , B is the amplitude of the meanders which varies 
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periodically in time as 



B = Bo + ecos{ut + (p) (14) 

The system can be fully mixing, i.e. any portion of the domain is definitely 
visited from any initial condition, in a certain portion of the parameter space 
{e,uj) where in particular one has cross stream transport. This system has 
two separatrices, with a spatial periodic structure, see Fig. [1], coinciding 
with the borderlines of the current. At very small e the chaotic layer is 
restricted to a limited region around the separatrices, and no cross stream 
mixing occurs. In order to have large-scale chaotic mixing, i.e. particles 
jumping across the jet from a northern recirculation to a southern one, and 



vice- versa, one needs the overlap of the resonances |jT6|, when e and uj are 
larger than certain critical thresholds. 

In Figure we report the OW indicator Ao(x) as function of the initial 
position. Incidentally, the HK indicator A+(x shows no significant differences 
from Aq and it is not shown. It is important to observe that both the OW 
and HK criteria are not able to detect the existence or not of a dynamical 
barrier, i.e. the two figures are practically indistinguishable. 

The finite time Lyapunov exponent 7r(x) is shown in Figure ^. The 
computation is done according to (0) with initial separation R{t) = S for 
all the particle pairs, without the "waiting procedure" (i.e. = 0). This 
is because we want to mimic a realistic situation of data analysis in which 
arbitrarily small separations cannot be attained. 

We note that, although the indicator is able to detect the jet core (the 
low FTLE value filament inside the chaotic current), the transition between 
the confined chaos regime and the large scale mixing regime is not observed 
(compare Figure ^ and Figure ^). The asymmetry of the FTLE map is 
due to the dependence of this indicator on the initial phase of the periodic 
flow. 

Let us consider the local FSLE as Ar(x), computed on the same trajec- 
tories and with the same initial separation of Figure ^. Figure ^ contains 
the results of the A^ maps, before and after the overlap of the resonances, 
in what we can call the Melnikov [|l7l and the Chirikov [jl6l regimes, respec- 
tively. The amplification factor r and the lower threshold 6 are chosen such 
that the upper threshold r ■ 5 is of the order of the jet width. In one case, 
black regions of zero FSLE values, i.e. particle pairs which never reach the 
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upper limiting separation, are located in the jetcore and in the centers of the 
recirculation: when chaos is still confined in the vicinity of the separatrices, 
no cross-stream transport is allowed. In the other case, after the separa- 
trix breaking has occurred, no zero FSLE values are present, indicating that 
particles can spread apart over any distance from any initial condition. 

Let us conclude this section by remarking that a Lagrangian diagnostic 
based on the FSLE shows is major skill in detecting dynamical barriers in 
the fiow. 



3.2 Geophysical data 

We now consider a geophysical example regarding Lagrangian motion on an 
isoentropic layer (i.e. at constant potential temperature) at a low strato- 
spheric level, in presence of the winterly polar vortex ||T^, characterized by 
a quasi-zonal robust jet stream. 

The typical fiow pattern is usually represented by means of stereographic 
maps of isoentropic Potential Vorticity (PV), see Figure ||, and the modulus 
of its gradient shown in Figure ^ In winter, the stratospheric PV can be 
considered as a quasi-conserved quantity over a timescale of about 2-3 weeks. 
It is also widely accepted in literature that the outer border of the polar 
vortex, usually identified by the maximum horizontal gradient of isoentropic 
PV 1 2^, can act as a strong barrier to meridional cross stream transport. 



The kinematics is remarkably similar to that of the previously discussed 
meandering jet model, if we imagine the latter as closed on itself in a cir- 
cular geometry. The Lagrangian data set we have used for the analysis 
account of about 10'* trajectory pairs, initially uniformly distributed over the 
southern hemisphere on the 475K isoentropic surface (lower stratosphere). 
Trajectories are calculated by means of the University of L'Aquila Trajec- 
tory Model 1^ using analyzed wind, pressure, and temperature fields 



from the U.K. Meteorological Office (UKMO) provided by the British 
Atmospheric Data Centre (BADC). Latitude coverage goes from poles to 
about tropics. The trajectories run from June 30th, 1997 up to a maximum 
observation time of 20 days. The initial distance between pair particles is 
Ri ~ 10 Km. 

Being the trajectory evolution simulated with wind fields relative to the 
southern hemisphere winterly season, we are observing a situation of stable 
polar vortex regime . 
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The local properties of the particle relative dispersion are obtained by 
computing the Okubo- Weiss eigenvalue Aq and the local FSLE A(r). Spatial 
derivatives of the velocity fields used for the calculations are estimated as 
finite differences over spatial grid steps of the order of 100 Km. 

The results are presented in the two-dimensional maps in Figures H and 
^. A further analysis (not shown) demonstrates that Aq does not change 
substantially in time and, as a consequence, it reproduces the essential fea- 
tures of the HK exponent A_|_. Neither of the exponents, as in the previous 
example, is able to detect the presence of a barrier to transport. 

On the other hand, the FSLE map (see Figure ||) detects the dynamical 
barrier as the region of vanishing FSLE values. The location of the barrier 
is in good agreement with the definition of polar vortex border based on 
geophysical considerations, e.g. the potential vorticity gradient shown in 
Figure 



4 Conclusions 

We have discussed several techniques proposed for describing dispersion in 
two-dimensional flows. In particular our analysis has been focused on the 
capability of these techniques to detect the presence of barriers to trans- 
port. By means of two examples of geophysical relevance, we have shown 
that Eulerian-based techniques, such as the Okubo- Weiss criterium and its 
generalization proposed by Hua and Klein, are not sensible to the presence 
of barriers. The Lagrangian finite-time Lyapunov exponent is, in principle, 
useful for describing space variations of the chaotic properties, e.g. in a 
quasi-integrable Hamiltonian system, but it is limited to small-scale proper- 
ties of dispersion. A recent non-linear Lagrangian diagnostics, based on the 
Finite Scale Lyapunov Exponent, is found to give the correct description of 
the presence of large-scale barriers. As final remark, we notice that the OW 
criterion has been recently shown to give poor information also in the case 
of fully developed turbulence |TB[: the probability distribution function of 



Ao, -P(Ao), for a typical 2D turbulent field is not sensitive to the presence of 
coherent structures, i.e. -P(Ao) is the same as for a Gaussian field. 

From a general point of view, it is not a surprise that purely Eulerian 
statistics, such as the OW quantity, are unable to predict the behavior of 
Lagrangian tracers. The presence of dynamical barriers is a fundamental 
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information about the transport properties of the flow and thus can be con- 
sidered as a good discriminatory for the diagnostics. It would be interesting 
to check the performance of the proposed methods on other geophysical flows. 

We thank the British Atmospheric Data Center for the UKMO data. We 
are grateful to B. Joseph, D. ludicone, B. Legras, R. Santoleri and G. Visconti 
for useful discussions. This work is partially supported by MURST (contract 
N. 9908264583). 
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Figure 1: Map of the stream function of the meandering jet model. The 
isohnes drawn in black represent the borderlines of the current. The jet 
core becomes a permeable barrier to cross-stream motion depending on the 
parameter values {u, e) of the time periodic perturbation. 
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Figure 2: Critical curve in the parameter space (a;, e) separating between 
"Melnikov regime" (below the curve), in which chaos is confined close to the 
separatrices, and "Chirikov regime" (above the curve), in which large-scale 
chaotic mixing occurs. Pi, (a;,e) = (0.1,0.3), and P2, (a;,e) = (0.4,0.3), are 
the two points in the parameter space discussed in the Lagrangian simula- 
tions, u! and e are adimensionalized with respect to ujq = 0.25 (pulsation of 
recirculating orbits next to the separatrices) and Bq — 1.2 (mean meander 
amplitude) . 
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a) 




Figure 3: Okubo- Weiss parameter Aq^^or the meandering jet system at the 
two parameter points Pi (a) and P2 (b) of Figure |^. 
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Figure 4: Finite time Lyapunov expor^t 7t(x) for the meandering jet system 
at the two parameter points Pi (a) and P2 (b). The number of particle pairs 
is 10000 with initial separation 5/L = 1.9 x 10~^ uniformally distributed on 
the periodic domain with spatial length L = 2-n jk (with k = 0.84). The time 
delay is r = tv/uq. 



a) 




2 4 6 

X 





0.0 0. 


5 - .0 


1 .5 2 




FSLE 






b) 






H il 






0.0 0.5 1 


.0 


.5 2 



FSLE 



Figure 5: The local FSLE Ar(x) for2Uie meandering jet system at the pa- 
rameter space points Pi (a) and P2 (b). The Lagrangian trajectories are the 
same of Figure 0. The amplification factor is r = 100. Only the particle 
pairs that reach a relative separation of r ■ 5 ~ lO^^L give a positive signal 
in terms of A(r). 
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Figure 8: Okubo-Weiss indicator Ao(x) for the analyzed wind fields corre- 
sponding to Figure computed at the seventh day of the Lagrangian simu- 
lation, July 6th 1997. 
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Figure 9: The local FSLE Ar(x) map for the isoentropic trajectory data set 
relative to the 475K layer. The initial distance between Lagrangian tracers 
is 5 ~ 10 km and the amplification factor is r = 10. 
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